Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon.

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

Genomic and pangenomic analyses provide insights into the population history and genomic diversification of bottle gourd.

Bottle gourd (Lagenaria siceraria (Mol.) Strandl.) is an economically important vegetable crop and one of the earliest domesticated crops. However, the population history and genomic diversification of bottle gourd have not been extensively studied. We generated a comprehensive bottle gourd genome variation map from genome sequences of 197 world-wide representative accessions, which enables a genome-wide association study for identifying genomic loci associated with resistance to zucchini yellow mosaic virus, and constructed a bottle gourd pangenome that harbors 1534 protein-coding genes absent in the reference genome. Demographic analyses uncover that domesticated bottle gourd originated in Southern Africa c. 12 000 yr ago, and subsequently radiated to the New World via the Atlantic drift and to Eurasia through the efforts of early farmers in the initial Holocene. The identified highly differentiated genomic regions among different bottle gourd populations harbor many genes contributing to their local adaptations such as those related to disease resistance and stress tolerance. Presence/absence variation analysis of genes in the pangenome reveals numerous genes including those involved in abiotic/biotic stress responses that have been under selection during the world-wide expansion of bottle gourds. The bottle gourd variation map and pangenome provide valuable resources for future functional studies and genomics-assisted breeding.

Antibody array-based proteome approach reveals proteins involved in grape seed development.

Grape (Vitis vinifera) is one of the most widely cultivated fruits globally, primarily used for processing and fresh consumption. Seedless grapes are favored by consumers for their convenience, making the study of seedlessness a subject of great interest to scientists. To identify regulators involved in this process in grape, a monoclonal antibody (mAb)-array-based proteomics approach, which contains 21,120 mAbs, was employed for screening proteins/antigens differentially accumulated in grape during development. Differences in antigen signals were detected between seeded and seedless grapes revealing the differential accumulation of 2,587 proteins. After immunoblotting validation, 71 antigens were further immunoprecipitated and identified by mass spectrometry (MS). An in planta protein-protein interaction (PPI) network of those differentially accumulated proteins was established using mAb antibody by immunoprecipitation (IP)-MS, which reveals the alteration of pathways related to carbon metabolism and glycolysis. To validate our result, a seedless-related protein, DUF642 domain-containing protein (VvDUF642), which is functionally uncharacterized in grapes, was ectopically overexpressed in tomato (Solanum lycopersicum "MicroTom") and led to a reduction in seed production. PPI network indicated that VvDUF642 interacts with pectin acetylesterase (VvPAE) in grapes, which was validated by BiFC and Co-IP. As anticipated, overexpression of VvPAE substantially reduced seed production in tomato. Moreover, S. lycopersicum colourless non-ripening expression was altered in VvDUF642- and VvPAE-overexpressing plants. Taken together, we provided a high-throughput method for the identification of proteins involved in the seed formation process. Among those, VvDUF642 and VvPAE are potential targets for breeding seedless grapes and other important fruits in the future.

The jet-like chromatin structure defines active secondary metabolism in fungi.

Eukaryotic genomes are spatially organized within the nucleus in a nonrandom manner. However, fungal genome arrangement and its function in development and adaptation remain largely unexplored. Here, we show that the high-order chromosome structure of Fusarium graminearum is sculpted by both H3K27me3 modification and ancient genome rearrangements. Active secondary metabolic gene clusters form a structure resembling chromatin jets. We demonstrate that these jet-like domains, which can propagate symmetrically for 54 kb, are prevalent in the genome and correlate with active gene transcription and histone acetylation. Deletion of GCN5, which encodes a core and functionally conserved histone acetyltransferase, blocks the formation of the domains. Insertion of an exogenous gene within the jet-like domain significantly augments its transcription. These findings uncover an interesting link between alterations in chromatin structure and the activation of fungal secondary metabolism, which could be a general mechanism for fungi to rapidly respond to environmental cues, and highlight the utility of leveraging three-dimensional genome organization in improving gene transcription in eukaryotes.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification.

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Identification of QTL associated with resistance to Phytophthora fruit rot in cucumber (Cucumis sativus L.).

Phytophthora fruit rot (PFR) caused by the soilborne oomycete pathogen, Phytophthora capsici, can cause severe yield loss in cucumber. With no resistant variety available, genetic resources are needed to develop resistant varieties. The goal of this work was to identify quantitative trait loci (QTL) associated with resistance to PFR using multiple genomic approaches and populations. Two types of resistances have been identified: age-related resistance (ARR) and young fruit resistance. ARR occurs at 12-16 days post pollination (dpp), coinciding with the end of exponential fruit growth. A major QTL for ARR was discovered on chromosome 3 and a candidate gene identified based on comparative transcriptomic analysis. Young fruit resistance, which is observed during the state of rapid fruit growth prior to commercial harvest, is a quantitative trait for which multiple QTL were identified. The largest effect QTL, qPFR5.1, located on chromosome 5 was fine mapped to a 1-Mb region. Genome-wide association studies (GWAS) and extreme-phenotype genome-wide association study (XP-GWAS) for young fruit resistance were also performed on a cucumber core collection representing > 96% of the genetic diversity of the USDA cucumber germplasm. Several SNPs overlapped with the QTL identified from QTL-seq analysis on biparental populations. In addition, novel SNPs associated with the resistance were identified from the germplasm. The resistant alleles were found mostly in accessions from India and South Asia, the center of diversity for cucumber. The results from this work can be applied to future disease resistance studies and marker-assisted selection in breeding programs.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato.

BACKGROUND: Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV). RESULTS: We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes. CONCLUSION: We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Chromosome-level phased genome assembly of "Antonovka" identified candidate apple scab-resistance genes highly homologous to HcrVf2 and HcrVf1 on linkage group 1.

Apple scab, a fungal disease caused by Venturia inaequalis, leads to losses in both yield and fruit quality of apples (Malus domestica Borkh.). Most commercial apple cultivars, including those containing the well-characterized Rvi6-scab-resistance locus on linkage group (LG) 1, are susceptible to scab. HcrVf2 and HcrVf1 are considered the main paralogs of the Rvi6 locus. The major apple scab-resistance loci Vhc1 in "Honeycrisp" and Rvi17 in "Antonovka," were identified in close proximity to HcrVf2. In this study, we used long-read sequencing and in silico gene sequence characterization to identify candidate resistance genes homologous to HcrVf2 and HcrVf1 in Honeycrisp and Antonovka. Previously published chromosome-scale phased assembly of Honeycrisp and a newly assembled phased genome of Antonovka 172670-B were used to identify HcrVf2 and HcrVf1 homologs spanning Vhc1 and Rvi17 loci. In combination with 8 available Malus assemblies, 43 and 46 DNA sequences highly homologous to HcrVf2 and HcrVf1, respectively, were identified on LG 1 and 6, with identity and coverage ranging between 87-95 and 81-95%, respectively. Among these homologs, 2 candidate genes in Antonovka and Honeycrisp haplome A are located in close physical proximity to the scab-resistance marker Ch-Vf1 on LG 1. They showed the highest identity and coverage (95%) of HcrVf2 and only minor changes in the protein motifs. They were identical by state between each other, but not with HcrVf2. This study offers novel genomic resources and insights into the Vhc1 and Rvi17 loci on LG 1 and identifies candidate genes for further resistance characterization.

Identification of genes associated with abiotic stress tolerance in sweetpotato using weighted gene co-expression network analysis.

Sweetpotato, Ipomoea batatas (L.), a key food security crop, is negatively impacted by heat, drought, and salinity stress. The orange-fleshed sweetpotato cultivar "Beauregard" was exposed to heat, salt, and drought treatments for 24 and 48 h to identify genes responding to each stress condition in leaves. Analysis revealed both common (35 up regulated, 259 down regulated genes in the three stress conditions) and unique sets of up regulated (1337 genes by drought, 516 genes by heat, and 97 genes by salt stress) and down regulated (2445 genes by drought, 678 genes by heat, and 204 genes by salt stress) differentially expressed genes (DEGs) suggesting common, yet stress-specific transcriptional responses to these three abiotic stressors. Gene Ontology analysis of down regulated DEGs common to both heat and salt stress revealed enrichment of terms associated with "cell population proliferation" suggestive of an impact on the cell cycle by the two stress conditions. To identify shared and unique gene co-expression networks under multiple abiotic stress conditions, weighted gene co-expression network analysis was performed using gene expression profiles from heat, salt, and drought stress treated 'Beauregard' leaves yielding 18 co-expression modules. One module was enriched for "response to water deprivation," "response to abscisic acid," and "nitrate transport" indicating synergetic crosstalk between nitrogen, water, and phytohormones with genes encoding osmotin, cell expansion, and cell wall modification proteins present as key hub genes in this drought-associated module. This research lays the groundwork for exploring to a further degree, mechanisms for abiotic stress tolerance in sweetpotato.

Genomes of cultivated and wild Capsicum species provide insights into pepper domestication and population differentiation.

Pepper (Capsicum spp.) is one of the earliest cultivated crops and includes five domesticated species, C. annuum var. annuum, C. chinense, C. frutescens, C. baccatum var. pendulum and C. pubescens. Here, we report a pepper graph pan-genome and a genome variation map of 500 accessions from the five domesticated Capsicum species and close wild relatives. We identify highly differentiated genomic regions among the domesticated peppers that underlie their natural variations in flowering time, characteristic flavors, and unique resistances to biotic and abiotic stresses. Domestication sweeps detected in C. annuum var. annuum and C. baccatum var. pendulum are mostly different, and the common domestication traits, including fruit size, shape and pungency, are achieved mainly through the selection of distinct genomic regions between these two cultivated species. Introgressions from C. baccatum into C. chinense and C. frutescens are detected, including those providing genetic sources for various biotic and abiotic stress tolerances.

Starch deficiency in tomato causes transcriptional reprogramming that modulates fruit development, metabolism, and stress responses.

Tomato (Solanum lycopersicum) fruit store carbon as starch during early development and mobilize it at the onset of ripening. Starch accumulation has been suggested to buffer fluctuations in carbon supply to the fruit under abiotic stress, and contribute to sugar levels in ripe fruit. However, the role of starch accumulation and metabolism during fruit development is still unclear. Here we show that the tomato mutant adpressa (adp) harbors a mutation in a gene encoding the small subunit of ADP-glucose pyrophosphorylase that abolishes starch synthesis. The disruption of starch biosynthesis causes major transcriptional and metabolic remodeling in adp fruit but only minor effects on fruit size and ripening. Changes in gene expression and metabolite profiles indicate that the lack of carbon flow into starch increases levels of soluble sugars during fruit growth, triggers a readjustment of central carbohydrate and lipid metabolism, and activates growth and stress protection pathways. Accordingly, adp fruits are remarkably resistant to blossom-end rot, a common physiological disorder induced by environmental stress. Our results provide insights into the effects of perturbations of carbohydrate metabolism on tomato fruit development, with potential implications for the enhancement of protective mechanisms against abiotic stress in fleshy fruit.

Transcriptomic and metabolic profiling of watermelon uncovers the role of salicylic acid and flavonoids in the resistance to cucumber green mottle mosaic virus.

Understanding the mechanisms underlying plant resistance to virus infections is crucial for viral disease management in agriculture. However, the defense mechanism of watermelon (Citrullus lanatus) against cucumber green mottle mosaic virus (CGMMV) infection remains largely unknown. In this study, we performed transcriptomic, metabolomic, and phytohormone analyses of a CGMMV susceptible watermelon cultivar 'Zhengkang No.2' ('ZK') and a CGMMV resistant wild watermelon accession PI 220778 (PI) to identify the key regulatory genes, metabolites, and phytohormones responsible for CGMMV resistance. We then tested several phytohormones and metabolites for their roles in watermelon CGMMV resistance via foliar application, followed by CGMMV inoculation. Several phenylpropanoid metabolism-associated genes and metabolites, especially those involved in the flavonoid biosynthesis pathway, were found to be significantly enriched in the CGMMV-infected PI plants compared with the CGMMV-infected 'ZK' plants. We also identified a gene encoding UDP-glycosyltransferase (UGT) that is involved in kaempferol-3-O-sophoroside biosynthesis and controls disease resistance, as well as plant height. Additionally, salicylic acid (SA) biogenesis increased in the CGMMV-infected 'ZK' plants, resulting in the activation of a downstream signaling cascade. SA levels in the tested watermelon plants correlated with that of total flavonoids, and SA pre-treatment up-regulated the expression of flavonoid biosynthesis genes, thus increasing the total flavonoid content. Furthermore, application of exogenous SA or flavonoids extracted from watermelon leaves suppressed CGMMV infection. In summary, our study demonstrates the role of SA-induced flavonoid biosynthesis in plant development and CGMMV resistance, which could be used to breed for CGMMV resistance in watermelon.

Genomic analyses reveal dead-end hybridization between two deeply divergent kiwifruit species rather than homoploid hybrid speciation.

Despite the importance of hybridization in evolution, the evolutionary consequence of homoploid hybridizations in plants remains poorly understood. Specially, homoploid hybridization events have been rarely documented due to a lack of genomic resources and methodological limitations. Actinidia zhejiangensis was suspected to have arisen from hybridization of Actinidia eriantha and Actinidia hemsleyana or Actinidia rufa. However, this species was very rare in nature and exhibited sympatric distribution with its potential parent species, which implied it might be a spontaneous hybrid of ongoing homoploid hybridization. Here, we illustrate the dead-end homoploid hybridization and genomic basis of isolating barriers between A. eriantha and A. hemsleyana through whole genome sequencing and population genomic analyses. Chromosome-scale genome assemblies of A. zhejiangensis and A. hemsleyana were generated. The chromosomes of A. zhejiangensis are confidently assigned to the two haplomes, and one of them originates from A. eriantha and the other originates from A. hemsleyana. Whole genome resequencing data reveal that A. zhejiangensis are mainly F1 hybrids of A. hemsleyana and A. eriantha and gene flow initiated about 0.98 million years ago, implying both strong genetic barriers and ongoing hybridization between these two deeply divergent kiwifruit species. Five inversions containing genes involved in pollen germination and pollen tube growth might account for the fertility breakdown of hybrids between A. hemsleyana and A. eriantha. Despite its distinct morphological traits and long recurrent hybrid origination, A. zhejiangensis does not initiate speciation. Collectively, our study provides new insights into homoploid hybridization in plants and provides genomic resources for evolutionary and functional genomic studies of kiwifruit.

Super-pangenome analyses highlight genomic diversity and structural variation across wild and cultivated tomato species.

Effective utilization of wild relatives is key to overcoming challenges in genetic improvement of cultivated tomato, which has a narrow genetic basis; however, current efforts to decipher high-quality genomes for tomato wild species are insufficient. Here, we report chromosome-scale tomato genomes from nine wild species and two cultivated accessions, representative of Solanum section Lycopersicon, the tomato clade. Together with two previously released genomes, we elucidate the phylogeny of Lycopersicon and construct a section-wide gene repertoire. We reveal the landscape of structural variants and provide entry to the genomic diversity among tomato wild relatives, enabling the discovery of a wild tomato gene with the potential to increase yields of modern cultivated tomatoes. Construction of a graph-based genome enables structural-variant-based genome-wide association studies, identifying numerous signals associated with tomato flavor-related traits and fruit metabolites. The tomato super-pangenome resources will expedite biological studies and breeding of this globally important crop.

SnRK1 kinase-mediated phosphorylation of transcription factor bZIP39 regulates sorbitol metabolism in apple.

Sorbitol is a major photosynthate produced in leaves and transported through the phloem of apple (Malus domestica) and other tree fruits in Rosaceae. Sorbitol stimulates its own metabolism, but the underlying molecular mechanism remains unknown. Here, we show that sucrose nonfermenting 1 (SNF1)-related protein kinase 1 (SnRK1) is involved in regulating the sorbitol-responsive expression of both SORBITOL DEHYDROGENASE 1 (SDH1) and ALDOSE-6-PHOSPHATE REDUCTASE (A6PR), encoding 2 key enzymes in sorbitol metabolism. SnRK1 expression is increased by feeding of exogenous sorbitol but decreased by sucrose. SnRK1 interacts with and phosphorylates the basic leucine zipper (bZIP) transcription factor bZIP39. bZIP39 binds to the promoters of both SDH1 and A6PR and activates their expression. Overexpression of SnRK1 in 'Royal Gala' apple increases its protein level and activity, upregulating transcript levels of both SDH1 and A6PR without altering the expression of bZIP39. Of all the sugars tested, sorbitol is the only 1 that stimulates SDH1 and A6PR expression, and this stimulation is blocked by RNA interference (RNAi)-induced repression of either SnRK1 or bZIP39. These findings reveal that sorbitol acts as a signal regulating its own metabolism via SnRK1-mediated phosphorylation of bZIP39, which integrates sorbitol signaling into the SnRK1-mediated sugar signaling network to modulate plant carbohydrate metabolism.

Characterization of the USDA Cucurbita pepo, C. moschata, and C. maxima germplasm collections.

The Cucurbita genus is home to a number of economically and culturally important species. We present the analysis of genotype data generated through genotyping-by-sequencing of the USDA germplasm collections of Cucurbita pepo, C. moschata, and C. maxima. These collections include a mixture of wild, landrace, and cultivated specimens from all over the world. Roughly 1,500 - 32,000 high-quality single nucleotide polymorphisms (SNPs) were called in each of the collections, which ranged in size from 314 to 829 accessions. Genomic analyses were conducted to characterize the diversity in each of the species. Analysis revealed extensive structure corresponding to a combination of geographical origin and morphotype/market class. Genome-wide associate studies (GWAS) were conducted using both historical and contemporary data. Signals were observed for several traits, but the strongest was for the bush (Bu) gene in C. pepo. Analysis of genomic heritability, together with population structure and GWAS results, was used to demonstrate a close alignment of seed size in C. pepo, maturity in C. moschata, and plant habit in C. maxima with genetic subgroups. These data represent a large, valuable collection of sequenced Cucurbita that can be used to direct the maintenance of genetic diversity, for developing breeding resources, and to help prioritize whole-genome re-sequencing.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis.

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C18 Δ9-elongase and a C20 Δ5-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ5-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Loss of OPT3 function decreases phloem copper levels and impairs crosstalk between copper and iron homeostasis and shoot-to-root signaling in Arabidopsis thaliana.

Copper (Cu) and iron (Fe) are essential micronutrients that are toxic when accumulating in excess in cells. Thus, their uptake by roots is tightly regulated. While plants sense and respond to local Cu availability, the systemic regulation of Cu uptake has not been documented in contrast to local and systemic control of Fe uptake. Fe abundance in the phloem has been suggested to act systemically, regulating the expression of Fe uptake genes in roots. Consistently, shoot-to-root Fe signaling is disrupted in Arabidopsis thaliana mutants lacking the phloem companion cell-localized Fe transporter, OLIGOPEPTIDE TRANSPORTER 3 (AtOPT3). We report that AtOPT3 also transports Cu in heterologous systems and contributes to its delivery from sources to sinks in planta. The opt3 mutant contained less Cu in the phloem, was sensitive to Cu deficiency and mounted a transcriptional Cu deficiency response in roots and young leaves. Feeding the opt3 mutant and Cu- or Fe-deficient wild-type seedlings with Cu or Fe via the phloem in leaves downregulated the expression of both Cu- and Fe-deficiency marker genes in roots. These data suggest the existence of shoot-to-root Cu signaling, highlight the complexity of Cu/Fe interactions, and the role of AtOPT3 in fine-tuning root transcriptional responses to the plant Cu and Fe needs.